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■ We study the existence of intraband and interband pairing in the extended two band and three-band 

Hubbard model. It is shown that including interband pairing significantly increases the supercon- 
ducting critical temperature in comparison with critical temperature of the intraband pairing. This 
increase allows for decrease of the nearest-neighbor interaction constant necessary for critical tem- 
perature in high temperature superconductors to the realistic values. 
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b ' 1. Introduction 
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^ High temperature superconductivity was discovered over fifteen years ago [1], but still today the mechanism of this 

phenomena is not uniquely established. The matter is made additionally more difficult by complicated band structure 
of the cuprates. The superconducting state is observed mainly in the Cu02 planes but it is created and influenced 
by the whole structure of cuprates, e.g. the apex oxygen [2-4] or other donor of localized level (Resonant Pinning 
Centers (RPC) [5]). Even description of the single superconducting Cu02 plane is rather complicated since in this 
plane the relatively lightly bound oxygen holes interact with strongly correlated copper holes. In the result we obtain 
hybridized oxygen and copper holes [6,7]. Interaction between the Cu02 planes and the interaction of copper and 
' oxygen holes within the Cu02 plane can be described by the two band model, in which strongly correlated levels are 
hybridizing with the broad conducting band [4,8-11]. The two band model can be used for superconductivity also for 
heavy fermion materials [12-14]. In most of the two-band models [4,8,11-13] it is assumed that the superconducting 
£^ , gap exist only in one hybridized band. Extension of this concept was presented in models in which there were two 
£f~) ' order parameters [10,15,16], but the order parameter in the band in which there is no Fermi level has very little 
, influence on the results and particularly on the critical temperature [10]. More interesting is the use of additional 
interband order parameters for superconductivity [17,18]. This concept is similar to the concept of interlayer order 
parameter used in the multilayer models [19-21]. In this paper we analyze the two-band model using both intraband 
and interband order parameters of superconductivity. We assume that the broad oxygen band is hybridized with the 
copper level what gives two bands separated with the hybridization gap. In a three-band model we analyze broad band 
interacting with two correlated copper levels ( 2>d x 2_ y 2 and 3d 3z 2_ r 2), what gives three hybridized bands. Including 
the 3d 3z 2_ r 2 level is stimulated by the experimental reports that there is a several percentage content of 3d 32 2_ r 2 holes 
in the electronic spectrum of layered cuprates (see [2,22]). In this paper we will concentrate mainly on the d-wave 
superconductivity. Series of experiments by ARPES [23] and SQUID [24] have shown that most of the cuprates have 
d-wave superconductivity. This kind of superconductivity is caused by the negative nearest-neighbor interaction W, 
, for which there is no unique explanation until now. One proposal was that in the three band model there is a significant 
difference; AV = V Z — V X , where V X (V Z ) is the nearest-neighbor interaction between copper orbital 3d x 2_ y 2 (3d 3z 2_ r 2) 
. i-h ! and the oxygen band 2p [25] . In the strong correlation limit the nearest-neighbor interaction W is equal to AV . 
^ ■ The problem is arising after we try to fit the value of negative AV to the critical temperature of the order of 100 K. 
This value for d-wave superconductor comes out to be unrealistically high [26] . Yet another source of the negative 
W can be the charge transfer instability [27,28]. Further on in this paper, to create the d-wave superconductivity, 
we assumed the existence of attractive nearest-neighbor interaction W, without going into details of its origin. In 
section 2 of this paper we propose the Hamiltonian describing interaction between the broad band (e.g. oxygen band 
2p) and the localized level (e.g. copper orbital 3d). This Hamiltonian is diagonalised by the standard procedure. 
To analyze superconductivity in hybridized bands we will use the Hartree-Fock factorization and the Green function 
formalism. In the result we obtain equations for the critical temperature of the d-wave and s-wave superconductor 
including both; intraband and interband pairing. In section 3 we analyze the dependence of critical temperature for 
d-wave superconductivity and nearest-neighbor interaction constant W, on the energy difference between center of 
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the broad band and the localized level. In addition we present in this section the influence of the interband pairing 
on values of critical temperature and nearest-neighbor interaction W in the three-band model. 



2. The Model Hamiltonian 

The Hamiltonian of the two-band model has the following form 

H = Ho + H T (1) 

where the kinetic part H is given by 
Ho = {s p — n) ^ 

J2 fe« + hx -) + ( £d - /*) 2 Ud " + f p d H ( d trP^ + h - c ) ( 2 ) 

icr ijcr icr ijcr 

with pf (pia) creating (annihilating) hole in the broad band, df (diu) creating (annihilating) hole in the localized 
level, n Pia — pf a Pi a is the hole number operator in the broad band, rid ia — d^dur is the hole number operator on the 
localized level, t pp is the broad band hopping integral, t pd is the integral of hopping between the broad band and the 
narrow level, e d is the energy of the narrow level, e p is the mean energy of the broad band. 
The potential part of the Hamiltonian Hi is given by 

Hi = U p J2 n Pj<r n Pj _ a +U d J2 n dia n di _ a + W ^ n d%a n P]a , , (3) 

jcr icr <ij>a,<T t 

where U p is the Coulomb repulsion in the broad band, U d is the Coulomb repulsion on the narrow level, and W is the 
attractive nearest-neighbor interaction between holes in the broad band and in the narrow level. 

We assume strong Coulomb correlation on the localized level (U d — > oo), which is taken as the zero energy level 
(e d = ) . This will change the kinetic part of the Hamiltonian to the following form 

H = (e p n Via + (Pt*Pi° + h - c -) + ( d t*Pi° + h.c.) - fij^ n ^ . ( 4 ) 

ia ija <ij><7 icr 

where the normalized hopping integral has the form 

tpd = tpdZ 1 / 2 , (5) 

with factor Z responsible for dynamic (hopping) electron correlation on the localized level. This factor can have one 
of the following forms 

Z = 1 — n d (mean-field approximation) (6) 

or 

Z = — (according to Gutzwiller approximation [29]). (7) 

1 — — d. 

1 2 

The kinetic part of Hamiltonian in the momentum space has the form 

H o = ^Z ( £fe ~ l l )PkaPk<? + (t P d,kdl a Pka + h.c.) - \i ^2 d~l a d kc r, (8) 

kcr kcr kcr 

where the dispersion relation for the broad band as can be seen from the geometry of the Cu02 plane, is given by 

£k = e P - 2t pp [cos(k x + k y ) + cos(fc x - k y )\ . (9) 
The interband hopping integral in the momentum representation has the form 

tpd.k = tpd ^ e lkS , (10) 
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where 5 is the distance between broad band atom (oxygen) and the nearest narrow level atom (copper). 

In the strong correlation limit (Ud — > oo) the second term in Hamiltonian (3) will disappear since there is no 
double occupancy on the copper atoms. The strong correlation, Ud, will contribute to the negative nearest-neighbor 
interaction W = V z — V x < 0, when we consider that in reality there are two orbitals on copper (3d x 2_ y 2 and 3e?3 Z 2_ r 2) 
with different strength of interaction V x , V z with adjacent oxygen orbitals 2p x , y (see [25]). In the result the potential 
part of the Hamiltonian (3) in the momentum representation will take on the form 

Hi = U p ^ pt a pt k _ a P-k>-*Pk>* + W kk ,d+ a pt k _ a P-k>-<jd kl(J i (11) 

kk'a kk'cr 

where 

W kk , = wY,z l{k ~ k ' )S = W{ lklk ,+ W ), (12) 

<5 



7fe = cos k x + cos k 



r\k 



cos k r — cos k„. 



(13) 



In pairing potential we neglected the triplet pairing terms. 

In this paper we do not analyze in details the source of the negative interaction W. We will only mention in 
here that in addition to the Weber's mechanism there is also a charge transfer from copper atoms to apex oxygen 
atoms laying directly above or beneath it. This charge transfer will create opposite charges on copper atoms and in 
plane oxygen atoms laying next to them, and result in the negative electrostatic W . This negative nearest-neighbor 
interaction will contribute to the d-wave superconductivity in addition to mentioned above Weber's mechanism. 

To analyze our two-band model first we have to diagonalize the kinetic part of our Hamiltonian given by Eq. (8) 
using the transformation between original (p, d) and new (a, (3) operators described by 



pi 
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where the dispersion relation for hybridized quasiparticles is given by 

1 



k 



(14) 



El 



The kinetic energy of the hybridized state has the form 



#o = ]T (25+ - »H a a ka + (E k - »Wjk„- 



(15) 



(16) 



ka 



ka 



In the next step we transform interaction Hamiltonian (11) using the same transformation of Eq. (14), which 
diagonalized the kinetic part of the Hamiltonian. To simplify further analysis we will neglect the weak Coulomb 
correlation in the broad band (U p = 0), because this repulsion in mean-field approximation does not influence the 
d-wave superconductivity [30, 31]. In the result we obtain 



H T = E W kkl L+a+ a +^-/3+J U± k a\_ a + uZ k P± ki _ a ) 

f uZ k ,(3-k'-a) (w + /afe',<7 +tJ k >Pk',a) 



kk' a 



iCt-k' 



(17) 



where 



u,„ = 



+ t 



and w+ = 



tpd.k 



pd,k 



(18) 



pd.k 



In Hamiltonianu (17) we have 16 combinations of operators a and 0. For each of them we use the Hartree-Fock 
factorization 



^ a k,a a -k,-a a -k',-<ra k ' 



: < a t,a a -k, 



> a-k' ,-aOtk' ,<r 



(19) 
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In the result, after omitting the constant terms (the last one in Eq. (19)), we obtain from Hamiltonian (17) the 
following equation 



where the functions fkd arc defined as 



/ a k,<7 \ 


-t- 






ftk,cr 









./i.,A, / M A* 

/fe.lAfe fkA^-k 

/ M A£ / M A* 

/fe.sAfe / fe , 2 A fe 



l -k.-c 



Pk,<7 

V #V„ J 



fk,i = u£u£,fk, 2 =^ k u k , f k<3 = UJ+U k , f kA = UJ k u+. 
The superconducting order parameter A^ is expressed as 

Afc = Ak,i + Afe_2 + Afe^ + Afc^, 

where 

Afe.i = ^ Wfefe'/fcM < alt^aifc/,-,, >, 



(20) 



(21) 
(22) 
(23) 



fc' 

A fc , 3 = J] Wfcfc' A',3 < <4 >(7 #±k',- ff >, 
fe' 

A M - J! ^fc' A'.* < Pk>,<, a -u,-o >■ 



(24) 



(25) 



(26) 



This order parameter is the sum of intraband parameters (A^^ and Afc.2), and interband parameters (A/^3 and 
Afc i4 ). Using Hamiltonian given by the sum of Eqs (16) and (20) in the equations of motion for Green functions 



e{{A;B)) £ = ([A,B} + ) + (([A,H]_-B)) £ , 



we obtain the following results 



S-E++/J, fk,iA k fkA^k 

f k ,iA* k S + / M A* 

/mA/c £-E k /fc. 2 A fe 

fkA^t / fe , 2 A* e + E^-v 



G(k,e) 



10 
10 
10 
1 



(27) 



(28) 



The solution of Eqs (28) is the Green function given by the matrix 

G(k,e) = 



z-K+V /fc,iA fc 
fk,iK e + E+-^ / M A* 



fkA^k 





where for example 





fkA^t 



G 12 (k,e) 



G 34 (k,e) 



f ki3 A k e-E k +^ f k , 2 A k 
f k , 2 A* k s + E k -fi 



/ M A fc [s 2 + (E k - 
D(e) 

f ka Ak [-e 2 + (El - M ) 2 
D(e) 



(29) 



(30) 



(31) 
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G 14 (k,e) 



J M A fc [-( £ - E~ + p)(e + g+ - M )] 



(32) 



G 32 (fc,e) — 

Denominator, D(e), for all these functions is given by 

D(s) = [f kA Al - £ 2 + (25+ - m) 2 ] [/£ 2 A 2 fe - £ 2 + (25 fe " - ^ 

+ [/kV* - (e - Sfe + + £ fe + - A 1 )] [/ fc 2 >4 Ag - (s + 25 fe " - M )( £ - £ fe + + M ) 
-2/ fe ,i/ fe , 2 / fe ,3/ fe ,4At - [e 2 - (£+ - M) 2 ] F - (25* - M) 2 ] 

Replacing in formulae (23)- (26) the mean values by the Green functions we obtain 

A fc ,i = -yV fcfc ,/ fc M / — hnG 12 {e,k')de, 



(33) 



(34) 



(35) 



A M = -Vw fe fc'/fc',2 / ^ImG 3 4(£,fc')&, 
A fc , 3 = -£W fcfc ,/fc',3 1 ^ImG 32 (£,fc')&, 



(36) 



(37) 



/(e) 

i / ImGi4(e, k')de. 



The total superconducting gap A^ is given by 



A k = -J2w kk , [ Mi m [ fk lGl2 (e,k') + f k ,, 2 G 34 (s,k') +f k ', 3 G 32 (s,k') + f k , A G 14 (e,k')]de. 



(38) 



(39) 



The symmetry of the gap A k depends on the symmetry of the pairing potential W kk >. From Eq. (12) one can see 
that the order parameter has to be written as 



A k = A s -f k + A d i] k . 



(40) 



Significant advantage of the current method lies in the fact that the superconducting state has pure d-wave and 
s-wave symmetry. In the two bands models with only intraband ordering [11-13] symmetry of the order parameter 
depended also on the hybridization factor (e.g. g k in [12,13]). The present definition of the order parameter (Eq. 
(40)) has the standard form for the s+d wave superconductivity. 

Inserting to Eq. (39) both the superconducting gap (40) and the effective potential (12), we obtain two separate 
equations describing two types of superconductivity 

- for the s-wave superconductivity 

1 = -W^ll I ^-Im[f k , A G* 12 (e,k') + f k ,, 2 G* 34 (e,k') +f k >, 3 G* 32 (e, k') + f k ,, 4 G* 14 (s, k')] de, (41) 

- for the d-wave superconductivity 

1 = -Wj2vlf ^-^[fk',iG* 12 (e,k') + f k , t2 G* 34 (s,k') +f k ', 3 G* 32 (E,k') + f k , A G* 14 (£,k')]ds, (42) 

k' 

where 
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G*Ae,k) 



Gij(e, k) 



(43) 



The Fermi level used in Eqs (41) and (42) is defined by the condition for carrier concentration in the narrow band 
"<* = E/ ^ImfK^Gu^fcO+w+Wfc [G 13 ( e ,fc) + G 31 ( e ,fe)] + (^) 2 G 33 (e,k)}de. (44) 

k 

Analytical solution of Eqs (41) and (42) for the nonzero gap (A& ^ 0, T < Tc) is very complicated. In the critical 
temperature; T — > Tc, — > 0, and the solution of Eqs (41) and (42) simplifies to the following forms 
- for the s-wave superconductivity 



k> 

f 2 , + f 2 , 

J k' ,3~ J k' ,4 



f^y tanh (I - M)) + tanh (f (££ - M )) 

anh (§(25+ - M )) + tanh (§(££ - M ))) 



for the d-wave superconductivity 

f 2 , 

J k' ,1 



(45) 



1 



fe' 



tanh 



J k' ,2 



2(b:,-m) 



tanh 



(f(^-M)) 



+ 



r, + f 2 , 

2(^+^-2^) 



(tanh (§(£+ - „)) + tanh (f (££ - M ))) 



(46) 



The above two equations together with Eq. (44) at Afc — » arc allowing for calculating the dependence of critical 
temperature on carrier concentration, which will be analyzed in the next section. 



3. Numerical results and discussion 

To compare Eq. (46) with analogous equation for the two-band model of superconductivity with order gap only in 
one band [11-13] we insert limiting values of fks = fk,3 = fk,4 = and 7^ into Eq. (46), obtaining the following 
equation 

where f k ,2 = E^t pd , k ' / [{E^f + i 2 pdk , . 

Including terms with fk t i, fk,3 and fk,4 in Eq. (46) allows for obtaining the critical temperature of 100 K at smaller 
values of nearest-neighbor interaction. Equation (47) gives this temperature at s = 2t p d, t pp = 0.15t p d, t p d = leV 
and rid = 0.75 for the nearest-neighbor interaction W = —0.7t p d, while from Eq. (46) the same temperature at the 
same band structure parameters can be reached already at W — —0.6t p d. This gives 15% reduction of the necessary 
interaction W. The main contribution to the reduction of W comes from the term with ft a, which represents the 
interband pairing (Eq. (20)). The other two terms (with f/-,i and /fc )3 ) have only small influence on W, e.g. the 
fk,i term, which is responsible for pairing in the upper sub-band (see Eq. (20)) decreases the value of the nearest- 
neighbor interaction W for only 1%. For this reason the upper sub-band was neglected in the previous papers with 
only intraband order parameter left [11-13]. For simplicity we assume that the interband hopping integral does not 
depend on the wave vector k (t p d.k = tpd)- 
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0.3 0.4 0.5 0.6 0.7 0.8 0.9 1. 

FIG. 1. Critical temperature Tc for the d-wave pairing as a function of the carrier concentration in the narrow band rid 
estimated from Eq. (46) (solid line) and Eq. (47) (dashed line), at e = 2t p d, t pp — 0.15t p d, t p a = leV, W = — 0.6t p d. 



At the same value of interaction W the model with interband pairing (Eq. (46)) will give larger values of the 
critical temperature than the model with only intraband pairing (Eq. (47)). Fig. 1 presents the critical temperature 
dependence on carrier concentration of copper estimated from Eqs (46) and (47). Both these curves were calculated 
for the same set of parameters; e = 2t p d, t pp = 0.15t p d, t p d = leV and W — — 0.679i p( j- One can see that the maximum 
of critical temperature is rising for about 50% when pairing takes place not only in one band. 

250 
200 
,_, 150 
E-h° 100 
50 


0.5 1 1.5 2 2.5 3 

s[eV] 

FIG. 2. Critical temperature Tc for the d-wave pairing as a function of energy difference between center of the broad band 
and the narrow level e with interband pairing (solid line) and without it (dashed line), at t pp = 0.15t p[ j, t P d = leV, W — —0.6t p d 
and rid = 0.75 



The influence of interband pairing on maximum of critical temperature and value of necessary nearest-neighbor 
interaction grows with decreasing energy difference between the localized level and the center of the broad band. Fig. 
2 shows dependence of the critical temperature on this energy difference e, with interband pairing (solid line) and 
without it (dashed line), at t pp = 0.15t pt j, t p d = leV, W = —0.6t p d and rid = 0.75. Dependence of the nearest-neighbor 
interaction necessary for reaching Tc — 100K on the energy difference e is given in Fig. 3 for t pp = Q.lbt p d, tpd = leV, 
rid = 0.75. Decreasing energy gap of hybridized bands; E£ — with decreasing e causes increase of the critical 
temperature and decrease of interaction W. 
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FIG. 3. Nearest-neighbor interaction W as a function of energy difference e with interband pairing (solid line) and without 
it (dashed line), at t pp — 0.15i p d, t p d = leV, n<j = 0.75and Tc = 100K. 



The above analysis will be even more important in the three band model describing the interaction between 
degenerated oxygen orbital O (2p) and two copper orbitals (3d x 2_ y 2 and 3(^2 _ r 2), see Gorski et. al. [26]. The 
kinetic part of this model is given by 



H = (Sp - M) E vX a Vka + (e x - A*) E d tka d xk* + {s z ~ fJ>) E d tka d ^a 
k,a k,a k,a 

E (rtxkd+kcrPka + k.C.) + ^zk d ^ ka Vkcr + h.C.) 



(48) 



k.a 



where 



. ni ( ■ k x a . k y a\ 2t pc i f . k x a , k y a\ ,..> 

t xk = 2t pd ^sm — - sin -J- j , t zk = — ^ ^sm — + sin -|- J , (49) 

d tka( d tka) ~ aie operators creating hole on the 3d x 2_ y 2(3d 3z 2_ r 2) level. 

The interaction part of Hamiltonian is given by Eq. (11) with d^ a {dka) being now the creation (annihilation) hole 
operators on the 2>d^, z 2_ r 2 level of copper. 

The effective pairing potential used by us previously (Eq. (14) in [26]) was only one from the 36 pairing terms coming 
from diagonalization of the model Hamiltonian. Although the chosen single term is dominating superconducting 
pairing, but as it was shown above, the other terms can be important as well, particularly terms with operators 
Xk, where there is a small energy difference between lower hybridized band pk and the middle band Xk- Equation 
describing order parameter and critical temperature of the d-wave superconductor in this the three-band model will 
have six terms and has the following form 

1 = -Wj2vl> [CfcM + Ck>,2 + C k ,, 3 + C k , tb + C k , t5 + C k , fi ], (50) 
fe' 

where coefficients C k >,i are given in Appendix. 

Numerical analysis of Eq. (50) shows that the nearest-neighbor interaction creating Tq rj 100K without interband 
pairing at e = 2t p d, t pp = 0, na = 0.81 is equal to W — — U.OAtpd. With the interband pairing, at the same structural 
parameters, this interaction drops to W = — 4.94£ pc ;- Replacing the oxygen level in the three band model by the oxygen 

band e k = e p — 2t pp cos( m t " ) — cos( "7 y ) with t pp = 0.5t p d (still e — 2t p d and rid = 0.75) reduces the interaction 

to W = —8.02t p d (see [26]). Including additionally interband pairing reduces this value further to W — —4.33t p d- 
Concluding; both effects of oxygen band broadening and interband pairing considered together allow for achieving 
small and realistic values of the nearest-neighbor charge-charge interaction. 



Acknowledgement 



We would like to express our gratitude to Prof. Karol I. Wysokinski from Marie Curie-Sklodowska University for 
suggesting the problem and for helpful discussions during the course of this work. 



Appendix 

Coefficients Ck>,% used in Eq. (50) are defined as follows 



Cfe' 



Cfe' 



(tpd 



25+ ) 2 tanh (f (25+ -/i)) 



{t pd E- k ,f tanh(f(i; fc --^)) 



Cfe' 3 — 



{t p dE k ,) 2 



(tanh (§(25+ - „)) + tanh (§(££ - M ))) 



3((2^)2 + ng / )((J#) 2 + ^) 



2(^+^-2/*) 



Cfc' 



fc',4 



(tanh (f(25+ - /x)) + tanh (f (25 fe " - M ))) 



3((^") 2 +^)((^) 2 + ^) 



2(^+^-2M) 



Cfe' 



(V^) 2 ( tanh (I ( E k> A*)) + tanh (f(25° - M ))) 



+ £fe°< - 



where 



Cfe/ fi — 



(t pd E k ,) 2 



(tanh (f(i5-, - „)) + tanh (f (25° - M ))) 



ng, ((^ fe -) 2 + ^) 



2(^fe" + ^ " 2/i) 



25+ = 1/2 



e± A /e 2 + 4fi2 



Sfe° = 



, , . k x a . k y a\ ^tpd ( . k x a , k y a 

= 4 sm ~z sm -77- H 7T- sm -tt" + sm -77- 

P d \ 2 2 / 3 V 2 2 
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